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Abstract. We present a new muon propagation Monte Carlo 
FORTRAN code MUM (MUons+Medium) which possesses 
some advantages over analogous codes presently in use. The 
most important features of the algorithm are described. Data 
on the test for algorithm accuracy are presented. Contribu- 
tions of different sources to the resulting error of simula- 
tion are considered. Selected results obtained with MUM 
are given and compared with ones from other codes. 



1 Introduction 

Propagation of muons in medium plays an important role for 
underground (-water, -ice) experiments with natural fluxes of 
high energy (HE) neutrinos and muons. Firstly, neutrinos are 
detected by muons which are born in i^N interactions and 
propagate a distance in medium from the point of interaction 
to a detector Secondly, muons which are produced in atmo- 
spheric showers generated by cosmic rays represent the prin- 
cipal background for neutrino signal and therefore their flux 
at large depths should be well known. Besides, the atmo- 
spheric muons deep under sea or earth surface are the only 
natural calibration source which allows to confirm correct- 
ness of the detector model by comparison experimental and 
expected detector response. For the muon propagation along 
with analytical methods one uses the Monte Carlo (MC) te- 
chnique which directly accounts for stochastic nature of the 
muon energy losses. There are several MC muon transport- 
ation algorithms currently in use (see, e.g. review in (Rhode 
and Carloganu, 1998)) but theoretical and experimental pro- 
gress makes to create new ones. 

Here we present a new MC muon propagation code MUM 
(MUonsH-Medium) written in FORTRAN. When working on 
MUM we aimed at creation of an algorithm which would: (a) 
account for the most recent corrections for the muon cross- 
sections; (b) be of adequate and known accuracy, i.e., does 
not contribute an additional systematic error which would 



exceed one from "insuperable" uncertainties (i.e. muon and 
neutrino spectra and cross-sections); (c) be flexible enough, 
i.e. could be easily optimized for concrete purpose to de- 
sirable and well understood equilibrium between CPU time 
and accuracy and easily extended for any medium and any 
correction for the cross-sections of the processes in which 
HE muon looses its energy; (d) be "transparent", i.e. pro- 
vide user with the whole set of data related to used models 
for the muon cross-sections, energy losses, etc.; (e) be as 
fast as possible. The MUM code has been developed for the 
Baikal experiment (Sokalski and Spiering, 1992; Balkanov 
et al., 2000) but we believe it to be useful also for other ex- 
periments with natural fluxes of HE muons and neutrinos. 



2 The basic features of the MUM algorithm ' 

To get finite CPU time Tcpu, the energy losses in any MC 
muon propagation algorithm have to be decomposed into two 
parts: muon interactions with fraction of energy lost v which 
exceeds some value Vcut, are simulated directly while the 
part of interaction with v < Vcut is treated by the approxi- 
mate concept of "continuum" energy loss. Setting Vcut too 
low one looses the speed (roughly, Tcpij oc v~^^) but set- 
ting it too high, one looses the accuracy. We did not fix 
Vcut in MUM, it may be set optionally within a range of 
10"'' < Vcut < 0.2, since the optimum value depends on 
the concrete case (Bugaev et al., 2000). 

An "absolute" energy transfer threshold AEcut in a range 
of 10 MeV < AEcut < 500 MeV can be used in MUM along 
with "relative" threshold Vcut to simulate the muon interac- 
tions within detector sensitive volume. It is important for 
deep underwater (-ice) Cherenkov neutrino telescopes (Sok- 
alski and Spiering, 1992; Balkanov et al., 2000; Barwick et 
al., 1991; Andres et al., 2000; Aslanides et al., 1999; Amram 
et al., 1999; Resvanis et al., 1994; Anassontzis et al., 2000), 
where the water or ice are used both as a shield which ab- 
sorbs atmospheric muons and as a detecting medium. 
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Detailed description can be found in (Sokalski et al., 2000) 
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Fig. 1. The simulated with MUM muon energy losses (markers) 
and model of energy losses as used by MUM (soUd hnes). Also 
results on the muon energy losses from (Lohmann et al., 1985) are 
presented (dashed lines). The plot corresponds to simulation in pure 
water with Ucut = 10~^. 



Fig. 2. Simulated with MUM distributions for the fraction of en- 
ergy lost in a single interaction v for muons with £ = 100 GeV 
(triangles) and E = \ PeV (circles) in comparison with correspond- 
ing differential cross-sections (Unes). The case for pure water and 
Vcut = 10~* is presented. The scales on Y-axis for E = 100 GeV 
and E = I PeV are different. 



The formulae for the cross-sections of muon interactions 
(bremsstrahlung, e+e~-pair production, photo-nuclear inter- 
action, knock-on electron production) are given in (Sokalski 
et al., 2000). The code does not use any preliminary com- 
puted files, all necessary data are prepared at the initiation 
on the base of several relatively short routines, which give 
cross-sections for the muon interactions and stopping-power 
formula for ionization. It allows user to correct or even en- 
tirely change the model for any type of the muon interaction. 

At each step we tried to avoid when possible any simpli- 
fications when computing/simulating the free path between 
two interactions, energy transfers, etc., or, at least, to track 
the error which comes from this or that kind of simplification. 
In most cases, to keep Tcpu at low level, roots for equations 
and values for functions are found at initiation, tabulated and 
then referenced when necessary by a interpolation algorithm 
which was carefully checked for each case to guarantee the 
high enough level of accuracy. 

Formally, MUM simulates the propagation of the muons 
with the energies up to 1 EeV but one should keep in mind 
that above 1 PeV uncertainties with muon cross-sections grow 
remarkably, some effects which expose at UHE (e.g., LPM 
effect) are not accounted in MUM. Three media are available 
instantly for the muon propagation with MUM, namely pure 
water, ice and standard rock. But any medium can be easily 
composed by user following examples which are given in the 
initiation routine. At its current version MUM represents an 
ID-algorithm which does not track the angular and lateral 
deviations of muons, but it is planned to be 3D-extended. 



3 The accuracy of the algorithm and an optimum set- 
ting of simulation parameters 

Fig.l shows results of an accuracy test which consisted of 
following. For a muon energy Eq the short distance Al was 
chosen and n muons were propagated through this distance 
with MUM 2. Then, (Eq - i X;r=i -^i)/^^ (where El 
is energy of i-th muon at the end of distance Al) represents 
energy losses simulated by the algorithm. In an ideal case 
they should be equal to ones which can be directly calculated 
by integration of differential cross-sections which represents 
the incoming model for given code but, since algorithm itself 
necessarily contributes an error which originates from appli- 
cation of numerical procedures which the code consists of, 
the simulated energy losses and incoming model for energy 
losses are not the same for any real code. The only case is 
presented in Fig.l but actually such test was performed both 
for water and standard rock with tJciu in a range from 10^'* 
to 0.05 (Sokalski et al., 2000). The difference between simu- 
lated energy losses and incoming model for energy losses for 
MUM does not exceed 1 % except for the case when v^ut > 
0.01 and ionization is treated as completely "continuous" 
process. It means that inner inaccuracy of MUM contributes 
to the resulting error much less than principal uncertainties 
with muons and neutrinos fluxes, cross-sections, etc. 

^The value / = n • A/ must be much greater than mean free path 
L{Eo) between two interactions with v > Vcut to obtain statisti- 
cally significant result but also must be small enough since there 
should be no stopped muons. 
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Fig. 3. Intensity of vertical atmospheric muon flux / at different 
deptfis D of pure water vs. Vcut as obtained by simulation with 
MUM. Muons were sampled according to sea-level spectrum from 
(Klimushin et al., 2000). Closed circles: knock-on electron produc- 
tion with fraction of energy lost v > Vcut is simulated; open circles: 
ionization is completely "continuous". Two horizontal solid lines 
on each plot show the flux intensity simulated with all muon cross- 
sections multiplied by a factor 1.01 (lower line) and 0.99 (upper 
line) for Vcut = 10~*. Dashed lines on plots for L) < 5 km corre- 
spond to intensity which was calculated for all energy loss treated 
as "continuous". Dash-dotted lines show the flux intensity simu- 
lated with muons sampled according to sea level spectrum (Gaisser, 
1990) and Vcut = 10"'*. Dotted Unes correspond to Vcut = 10~* and 
cross-section for absorption of a real photon at photo-nuclear inter- 
action parameterized according to (Breitweg et al., 1999) instead of 
parameterization (Bezrukov and Bugaev, 1980,1981) which is the 
basic in MUM. 



Fig.2 demonstrates the accuracy of simulation for fraction 
of energy lost v for different kinds of muon interaction. Sim- 
ulated distributions are plotted along with functions for dif- 
ferential cross-sections da/dv. Again, only small part of the 
data is shown in the plot but agreement between simulated 
distributions and predictions is not worse for other media and 
other muon energies which has been tested carefully. 

Results on simulation of atmospheric muon vertical flux 
at different depths in the pure water vs. Vcut are shown in 
Fig. 3. Simulations were performed for 2 atmospheric muon 
surface spectra; with knock-on electron production included 
in simulation of energy losses and treated as totally "contin- 
uous"; with 2 different parameterizations for photo-nuclear 
interaction; with all muon energy losses multiplied by 0.99, 
1.00 and 1.01 (which corresponds to the most optimistic ev- 
aluation for uncertainties with the muon cross-sections of 
1% (Rhode and Carloganu, 1998; Kokoulin and Petrukhin, 
1991,1999)). The general conclusion is as follows: the prin- 
cipal uncertainties when computing the atmospheric muon 



flux at large depth are ones for the muon cross-sections. In- 
fluence of set value for Vcut in a range 10^^-10^^ is much 
less and, in principal, in a "ideal muon propagation code" one 
could set Vcut = 10~^ which allows calculations to be rather 
fast without remarkable influence on the result. Also the 
ionization energy losses can be treated as completely "con- 
tinuous" (which saves Tcpu with a factor of ~2). But, as 
MUM's own accuracy (in the sense of reproducing the muon 
energy losses) becomes worse than 1% for Vcut > 10~^ 
if ionization is excluded out of simulation and for Vcut > 
5-10~^ if knock-on electrons are simulated (Sokalski et al., 
2000) we conservatively affirm Vcut = 5 -10"^ and knock-on 
electron production included in simulation as a optimum set- 
ting of parameters for simulation the atmospheric muon flux 
at large depths with MUM. With such setting the proportion 
of Tcpu which is necessary to get the same statistics with 
muon propagation algorithms MUM, PROPMU (Lipari and 
Stanev, 1991) and MUSIC (Antonioh et al., 1997) is approx- 
imately 1 : 10 : 600 (note tiiat MUM, in contrast botii to 
PROPMU and MUSIC, is ID algorithm). 

We did not investigate specially the influence of simulation 
parameters on the results for the muon flux originated from 
neutrino. Generally, intensity of the muon flux I^*^ which 
accompanies the neutrino flux in a medium is proportional 
to the muon range, and, consequently, 7^*^ oc {dE/dx)~^ 
(in contrast to atmospheric muons whose flux at the large 
depths depends more sharply upon muon energy losses - see, 
e.g.. Fig. 3). That means that an error for simulated flux of 
muons produced by neutrino is proportional to an error in 
muon energy losses. So, the setting of parameters described 
above fits even better for propagation of muons originated 
from neutrino. 



4 Comparison to other muon propagation algorithms 

Fig. 4 shows survival probabilities (fractions of muons which 
have survived after propagation of distance D) vs. distance 
of propagation in pure water as computed with MUM, PRO- 
PMU and MUSIC for a set of muon energies from 500 GeV 
to 30 PeV. There are no statistically significant differences 
between MUM and MUSIC but survival probabilities com- 
puted with PROPMU are noticeably higher at energies E < 
30 TeV (up to a factor of 6) and become less at > 30 TeV. 

In Fig. 5 differential spectra for vertical atmospheric mu- 
ons at different depths in pure water are presented as simu- 
lated with MUM, PROPMU and MUSIC. Muons at the sur- 
face were sampled according to spectrum (Klimushin et al., 
2000). Okada (Okada, 1994) and KBS (Klimushin et al., 
2000) parameterizations for deep underwater muon spectra 
are shown, as well. MUSIC and MUM give almost the same 
results because survival probabiUties for muons in pure wa- 
ter are the same when simulating with MUSIC and MUM. 
music's and MUM's spectra coincide with KBS paramete- 
rization which is based on the same sea-level muon spectrum 
as was used for simulation and on muon propagation with 
MUM. Okada parameterization is lower than KBS, MUM 
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Fig. 4. Survival probabilities vs. distance of propagation in pure 
water as simulated with MUM (solid lines), PROPMU (dashed 
lines) and MUSIC (circles). Figures near curves indicate initial en- 
ergy of mono-energetic muon beam as follows: 500 GeV (1), 1 TeV 
(2), 3 TeV (3), 10 TeV (4), 30 TeV (5), 100 TeV (6), 300 TeV (7), 1 
PeV (8), 3 PeV (9), 10 PeV (10), 30 PeV (11). 

and MUSIC results (up to 18% in terms of integral muon 
flux at = 1 km) at relatively shallow depths and becomes 
higher at Z) > 5 km because it is based on rather hard surface 
muon spectrum with index 7 =2.57 which leads to a deficit 
for low energy muons comparing to KBS parameterization. 
Simulation with PROPMU produces the muon spectra which 
i) are significantly higher (31%, 30%, 27% and 17% in terms 
of integral muon flux at the depths D = I km, 3 km, 6 km 
and 10 km, correspondingly) and ii) are expanded to the low 
energies. It is in qualitative agreement with results on survi- 
val probabiUties presented in Fig. 4. 

5 Conclusions 

We have presented the muon transportation algorithm MUM 
and have given selected results obtained with it in compari- 
son with ones obtained with analogous codes. We consider 
the current version of MUM as a basis for the further devel- 
opment. The code is available by request. 
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